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^1 '■ Abstract. On using the known equivalence between the presence of a position-dependent 

^< | mass (PDM) in the Schrodinger equation and a deformation of the canonical commutation 

I/"") ' relations, a method based on deformed shape invariance has recently been devised for genera- 

ting pairs of potential and PDM for which the Schrodinger equation is exactly solvable. This 
approach has provided the bound-state energy spectrum, as well as the ground-state and 
<**j \ the first few excited-state wavcf unctions. The general wavefunctions have however remained 

Oh 1 unknown in explicit form because for their determination one would need the solutions of 

^ (— i a rather tricky differential-difference equation. Here we show that solving this equation may 

be avoided by combining the deformed shape invariance technique with the point canonical 
transformation method in a novel way. It consists in employing our previous knowledge of 
the PDM problem energy spectrum to construct a constant-mass Schrodinger equation with 
similar characteristics and in deducing the PDM wavefunctions from the known constant- 
£N1 ' mass ones. Finally, the equivalence of the wavefunctions coming from both approaches is 

checked. 
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1 Introduction 
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There exists a wide variety of physical problems in which an effective mass depending on the 
$_i ' position is of utmost relevance, such as effective interactions in nuclear physics [I], carriers and 

impurities in crystals [2j, quantum dots [3], quantum liquids j4j, semiconductor heterostruc- 
tures [5], and physics in neutron stars [6|. Position-dependent masses (PDM) also hold out to 
deformation in the quantum canonical commutation relations or curvature of the underlying 
space [9]. Furthermore, they may also appear in Hermitian Hamiltonians equivalent to 
"PT-symmetric or pseudo-Hermitian ones |10|. [IT] , whose study is a topic of current considerable 
interest [T2] . 

It is worth stressing that in addition to their position dependence responsible for noncom- 
mutativity with the momentum [13], some PDM may present unusual characteristics, such as 
discontinuities, singularities or vanishing at one or both end points of the interval. These pro- 
perties may create some problems, which are not present in the constant-mass case, but are 
dealt with in the literature. For instance, a vanishing of the PDM, which often occurs in curved 
spaces [3 El [9] with applications to quantum dots [H] , imposes some supplementary condition 
on bound-state wavefunctions [15] . 

Much attention has recently been devoted to finding exact solutions to Schrodinger equations 
in a PDM context due to their usefulness in physical applications. Such studies use all kinds 
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of methods known for solving constant-mass Schrodinger equations or an extension of them. 
References on the subject can be traced, for example, through [IB]. In the present paper, we 
quote only the ones that are directly relevant to what we are going to deal with. 

On building on the previously noted equivalence between the presence of a PDM and a de- 
formation of the commutation relations [7J, we have recently devised a method for generating 
pairs of potential and mass for which the Schrodinger equation is exactly solvable [15] . For 
such a purpose, we have used an approach inspired by a branch of supersymmetric quantum 
mechanics |17[ \TE[ [19] , whose development dates back to that of quantum groups and g-algebras 
[20[ [2"T] [22] [23] [2^1 [25] 126] . It consists in considering those one-dimensional potentials that are 
translationally shape invariant for a constant mass and in deforming the corresponding shape 
invariance condition in such a way that it remains solvable. 

This type of method has easily provided us with exact results for the bound-state energy 
spectrum and for the ground-state wavefunction. However the derivation of the excited-state 
wavefunctions in explicit form has turned out to be far more tricky because they are expressed 
in terms of polynomials satisfying some differential-difference equation, whose solution is rather 
difficult in general form. Up to now we have only been able to directly solve it in the case of 
the <i-dimensional radial oscillator [27]. This result has then been extended to the Morse and 
the D-dimensional Coulomb potentials [IB] by using a point canonical transformation (PCT) 
analogous to that relating the constant-mass problems |28 [ I29 [ 130]. In doing so, we have taken 
advantage of the fact that the three potentials belong to the so-called Natanzon confluent 
potential class [31] G2] . 

The remaining point at issue is therefore the explicit form of the excited-state wavefunctions 
for those potentials of [15] that do not belong to this class. It is the purpose of the present 
paper to answer such a question. To this end, it will prove expedient to combine the deformed 
shape invariance (DSI) technique with the PCT method in a novel way. Since we already know 
the energy spectrum of each PDM problem, we may try to find a constant-mass Schrodinger 
equation giving rise to a similar type of spectrum and to devise a PCT mapping the former 
problem onto the latter. It will then be a simple task to employ the reciprocal PCT to derive 
the PDM Schrodinger equation wavefunctions from the knowledge of the constant-mass ones. 
As we plan to show, this combined method works for all pairs of potential and mass considered 
in [15] , including those already solved in |16[ [27] . 

This paper is organized as follows. In Section the application of the DSI method to 
PDM Schrodinger equations is reviewed and the previously obtained results summarized. In 
Section [3] the use of the PCT method is explained by means of two detailed examples and the 
wavefunctions are listed for all the potential and PDM pairs. Finally, Section O contains the 
conclusion. 



2 Deformed shape invariance method 

In one-dimensional nonrelativistic quantum mechanics, one may deform the conventional cano- 
nical commutation relation [x,p] = i, where p = —id/dx and H = 1, into [7J 

[x,tt] = if(a;x), vr = -iy/ f(a;x)-^y/ f(a;x). (2.1) 

Here we assume that the deforming function f(a;x) is real, positive, and depends on a set of 
real parameters a, in such a way that f(a;x) — > 1 in the a — > limit. In the following, it 
proves convenient to write f(a;x) = 1 + g(a;x). Note that from (|2.ip . it results that both x 
and 7r are Hermitian operators. 

By substituting it 2 for p 2 in the conventional Schrodinger equation 

Hil> n (x) = [p 2 + V{x)]^ n (x) = E n ^ n (x), (2.2) 
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we arrive at a deformed equation 

H^^\x) = [vr 2 + V{x)]^\x) = E^^\x). (2.3) 
The latter may be re-interpreted as a PDM one, 

1 + VM^x)) ^(x) = E^^(x), (2.4) 



dx M(a; x) dx 
with a mass and an effective potential given by 

M(a; x) = 1 V efi (a; x) = V{x) - ±/(a; x)f"(a; x) - ]f' 2 (a; x), (2.5) 

f z {cx;x) 2 4 

respectively [7J \15\ I16j . Here a prime stands for derivative with respect to x. In (|2.2p . we 
have taken units wherein the constant mass mo = 1/2, while in (|2.4j) . M(q;x) denotes the 
dimensionless part of the mass function m(a; x) = moM(at; x). 

It should be noted that the ordering chosen for the noncommuting momentum and mass 
operators in the PDM Schrodinger equation (|2,4p is that of BenDaniel and Duke [33], for which 
some physical arguments have been put forward [34]. Other choices maintaining Hermiticity of 
the kinetic energy operator may be taken care of by adopting the von Roos ansatz, depending 
on two independent ambiguity parameters [13]. In such a case, the effective potential, defined 
in (|2.5p . contains an additional contribution coming from those parameters (see, e.g., |35j for 
a discussion of this topic) . 

The DSI method [15] considers the Hamiltonian H^ a \ defined in equation (j2.3l) . as the first 
member Hq°^ = of a hierarchy of Hamiltonians 

i 

= A+(«,A 4 )A-(«,A 4 )+^e„ i = 0,1,2,..., 

where the first-order differential operators 

d 



A ± (a,X l ) = ^y/f(a;x) — y/f(a;x) + W(Xi;x) 

satisfy a DSI condition 

A-(a,\ i )A+(a,\ i ) = A+(cx,\ i+1 )A-(a,\ i+1 ) + e i+1 , i = 0,1,2,..., (2.6) 

and ej, i = 0, 1, 2, . . ., are some real constants. Solving equation f|2.6|) means finding a super- 
potential W(X;x) (depending on some set of real parameters A), a deforming function f(a;x) 
and some constants Aj, Cj, i = 0, 1, 2, . . ., with Ao = A, such that 

V(x) = W 2 {\- x) - f{a; x)W'{\; x) + e (2.7) 

and 

W 2 {X l ;x) + f(a;x)W'(\ i ;x) 

= W 2 (X i+1 ;x)- f(a;x)W'(X l+1 ;x) + e i+1 , i = 0,1,2,.... (2.8) 

On starting from the known superpotentials of translationally shape-invariant potentials [17] . 
we have been able in most cases to maintain the solvability of equations (|2.7p and (|2.8[) for / / 1 
by a procedure detailed in [15] . The resulting deforming functions and deformed superpotentials 
are listed in Tables [T] and [2] for the shifted harmonic oscillator (SHO), radial harmonic oscillator 
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Table 1. Potentials, deforming functions f(a;x) = 1 + g(a;x) and effective potentials. 
Type V(x) 9{ol\x) V eii (a; x) - V(x) 

SHO a \uj 2 {x-^f ax 2 + 2(3x -2(ax + p) 2 - 5 2 

-oo<x<oo a > j3 2 > 0, / -22* <5 = - (5 2 

SHO fe \oj 2 {x-^f ax 2 + 2(3x -2a 2 (x - f) 2 - 5 2 

-oo<x<oo a > f3 2 > 0, (3 = 5 = Ja (l - 



RHO Vx 2 + ^il ax 2 -2a 2 x 2 -a 

< x < oo a > 

C -f + ^+li ax 

< X < OO a > 

M B 2 e~ 2x - B{2A + l)e~ x ae~ x -\a 2 e~ 2x - \ae~ x 

— oo < x < oo, A, B > a>0 

E a A(A — 1) csch 2 x — 2B cothx ae~ x smhx ae~ 2x (5 — \ae~ 2x ) 

0<x<oo, A>|, #>A 2 -2<a/0 <5 = 1 + f 

E;, A(A — 1) csch 2 x — 2B coth x ae~ x sinh x — 3e~ 4x 

< x < oo, A > §, B > A 2 a = -2 

PT A(A — 1) sec 2 x asin 2 x — a — 2a(a — 1) sin 2 x 

-§<£<§, A>1 -Ka^O + 3a 2 sin 4 x 

S [A(A - 1) + B 2 ] sec 2 x asinx § sinx - ^(1 - 3sin 2 x) 

- B(2A - l)secxtanx < \a\ < 1 

-f < x < f, A- 1 > 5 > 



RM j4(v4 — 1) esc 2 x + 2.B cot x sinx(acosx + /3sinx) ( 1 + Ih (a sin 2x 
0<x<vr, A>1 VTTP > /3 > -1 -/3cos2x) 



| (a 2 - /3 2 ) cos4x 
|a/3 sin 4x 



+ Ua 2 + (3 



(RHO), radial Coulomb (C), Morse (M), Eckart (E), Poschl-Teller (PT), Scarf I (S) and Rosen- 
Morse I (RM) potentials, respectively. It is worth noting that the sets of parameters a. and A 
contain either one element (denoted by a and A) or two elements (denoted by A, fi and a, (3) 
and that the radial harmonic oscillator and Coulomb potentials actually act in a D-dimensional 
space, where the angular momentum quantum number I is related to the parameter L through 
the equation L = I + Furthermore, the Scarf II potential is missing from the lists because 
no positive function / on the whole real line has been found. In contrast, the Rosen-Morse II and 
generalized Poschl-Teller potentials give rise to some acceptable functions /, but are omitted 
from Tables [Q and [2] for a reason explained below. 
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Table 2. Superpotentials and corresponding parameters. 



Type W(\,x) A A 



SHO afe Xx + u A = \(a + A), A = Vw 2 + a 2 X n = X + na 



RHO 7 + A = —L — 1 A n = A- n 



a 2 fj, n = n + na 
| + A = — L — 1 X n = A - n 



A 


A = 


2(a + Aj, A 


I 1 = 


^ ~~ 2A 


\ - 

A — 


T 1 


I 1 = 


1 / „ 1 A \ A 

2 (a + A), A 


A = 


—L - 1 


V = 


2Z+a\ 
2X 


A = 


-i(a + A), Z 


\i = 


1 (B(2A+1) 

2 \ A 


A = 


-A 


V- = 




A = 


|(1 + a + A) 



\[i+2n/3\+n 2 a/3 
\+na 



2Z+a\(2n+l)-an 2 
2(A-n) 



M Ae- X + ^ A = -i(a + A), A = V4S 2 + a 2 A n = A - na 

\ _ 2\{n-n)+n 2 a 

) — 2(A-na) 

E aj 6 A coth x + u X = — A X n = X — n 

A/i— ^an(2A— n) 
^ n = A^n 

PT Atanx A=i(l + a + A) X n = X + n(l + a) 



A = y/(1 + a) 2 + AA{A - 1) 
A tan x + /x sec x A = ^(1 + A + + A_) A n = A + n 

/i = ^(a — A + + A_) u n = fi + na 



A -± = ^(1=F«) 3 + C7 ± (C ± -1) 
C± = A ± 5 

RM Acotx + zu A = — A A n = A — n 

B 1 _ A/i— |an(2A— n) 

^ — _ A - 2° /^n — X^n 

From the constants ej, the energy eigenvalues are determined through the relation 

n 

The ground-state wavefunction, annihilated by A~ (ex, A), can be written as 

4") (a;) = 4 a) (\; x) oc - 1 exp (- f W ,j^^ dJ\ (2.10) 



and the excited-state wavefunctions can in principle be obtained recursively by acting with 
A+(a,X), 

= Vffi (A;x) oc 4+(a, A)#(A i; x). (2.11) 

Equation (|2.9p only provides solutions to the bound-state energies of if the corresponding 
wavefunctions (|2.10p and (|2,lip are physically acceptable. As stressed in [To] , this imposes not 
only that they are square integrable on the (finite or infinite) interval of definition (x±,X2) 
of V(x), as for the conventional Schrodinger equation, but also that they ensure the Hermiticity 
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of (or that of tt). The latter condition translates into 

| ^4") (A; x)| 2 /(q; x) = —. = ' = >0 for x — > x\ and x — > X2, 

y/M(ct; x) 

which imposes an additional restriction whenever f(a;x) — > oo (or M{a;x) — ► 0) for x — > xi 
and/or x — ► X2- It is actually this condition that is not satisfied by the square-integrablc 
wavefunctions obtained for the Rosen-Morse II and generalized Poschl-Teller potentials and 
explains their absence from Tables [T] and [2j 

(at) 

The resulting bound-state energies E n , corresponding to the potentials and deforming func- 
tions of Table [H are listed in Table [3j This table illustrates the strong influence that the de- 
formation or mass parameters may have on the spectrum. This is particularly striking for the 
Coulomb potential, whose infinite number of bound states for a constant mass is converted into 
a finite one, and for the Eckart potential, for which a finite number of bound states becomes 
infinite in the case a = —2. We plan to come back to these features in Section [3] and to provide 
there a simple derivation of the bound-state number. 

Table 3. Bound-state energy spectra. 

Type E n a) 

SHO a , b (2n + 1)A + n 2 a + d 2 - ( W f + l + f a(3 ) 2 0,1,2,... 



X+na 

RHO -2Au - aX + u - 4(aA - u)n + 4cm 2 0,1,2,... 

p 1 / 2Z+aX(2n+l)-an 2 \ 2 n 1 

^ -4 I J u, i,...,n max 



<ax + |A|(2n max + l) < f 
< K iax + l) 2 + |A|(2n max + 3) 



M -1 ( 2A( T-ir a ) 0,l,...,n max 

(2|A| +n max a) < 2|A|u 

< (n max + 1)[2|A| + (n max + l)a] 

"(A - n) 2 - ( 2A %T-n)"" } ) 2 °' X > ■ • • > n — 

+ a[(2n + 1)A - n 2 ] (| A| + n max ) 2 < |A|(2 ^ +a|A|) 

< (|A| +n max + l) 

_( A _ n )2_( A M+ n(2A-n) X 2 0,1,2,... 



2d 

2 



Eh 



-2[(2n + l)A-n 2 ] 

PT (A + n) 2 -a(A-n 2 ) 0,1,2,... 

S a (A + n) 2 - a(2n + l)u - a 2 n 2 0,1,2,... 

(A-n) 2 -( %7Jg-") ) 2 0,1,2,... 
- /3(2n + 1)A + (3n 2 



RM 



a A misprint has been corrected in the Appendix of [15] . 

The corresponding wavefunctions are listed in Tabled! where they are given in terms of nth- 
degree polynomials P n (\;y) in some new variable y, such that Po(X;y) = 1. As a consequence 
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of the recursion relation (|2.1ip . those polynomials satisfy differential-difference equations, given 
in Table From such equations, it is clear that if the ground-state and first few excited-state 
wavefunctions are easy to obtain, the same is not true for the remaining ones. 



Table 4. Bound-state wavefunctions resulting from the DSI approach. 



Type 



V4°°(A;x) 



SHO a , fe 


/ 


RHO 


X 


C 


X 


M 


f 


E Q 


(< 


E 6 


(c 


PT 


(< 


S 


/ 


RM 


/ 



(A n + a) 

2a exp 



-^^arccot^)p n (A,/x;y) 



l A l/ 2. P n (A,/i;y) 

2A n -(2 Mn + l) a 

2a e ^ x P n (\,^y) 



(l+a)X n — IHi 1 An+Mn . . 

I 25 2(cothx-l) 25 P n (A,^i;y) cothx 



(coth x — 1 csch x exp 

A n A n +l+q 



/I 

coth x- 



T ) P n (A,//;y) 



An — aj^Tt 1 



A + M 



|(|A n |+l)/ ■ „\|A„| 



sinx)' ™l exp 



2/X n +a|A n | 

2<5 



arccot cotx s + 2 j p n (A, ^; y) 



cothx 
tanx 
sinx 
cot x 



1 + 



Table 5. Differential-difference equations satisfied by the polynomials. 



Type Equation for P n (A;y) 







RHO 


P n+ i(A,/i;y) = ( 




x P n (Ai,//i;y) 


C 


P n+ i(A,/i;y) = | 


M 


P n+ i(A,/i;y) = | 


E a ,6 


P n+ i(A,/i;y) = | 


PT 


Pn+i(A;y) = (- 


S 


P n+ i(A,//;y) = | 




x P n (Ai,/ii;y) 


RM 


P n+ i(A,/i;y) = | 



-(1 + 2/?y + ay 2 )^| + (A n+ i + A)y + + (i) P n (Ai, ^; y) 



-2y(l + ay) + A n+ i + A + n + (// n+1 + \i + na)y) 



y(y + a )-^ + + A)y + /x n+ i +//J P„(Ai,/ii;y) 

y(l + ay)^ + (A n+ i + A)y + // n+ i + fij P„(Ai,^i;y) 

(y 2 + ay - l - + (A n +i + A)y + /i n +l + Pn(Ai,/ii;y) 

1 + (1 + a)y 2 ] | + (A n+1 + A)y) P„(A i; y) 

-(1 + ay)(l - y 2 )^ - ny(l + ay) + (A n+ i + A)y + fi n+ i + fi 



(y 2 + ay + 1 + (3)± + (A n+ i + A)y + /x n+ i + fi\ P n (Xt,m;y) 



To obtain all the bound-state wavefunctions in explicit form, we shall proceed in Section 
to combine the DSI approach with the PCT method. 
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3 Point canonical transformation method 
3.1 General method 

As shown in [35], a constant-mass Schrodinger equation 
d 2 \ 

--j-^ + U(u) J <t> n {u) = s n <t) n {u), (3.1) 

for some potential U(u) defined on a finite or infinite interval, can be transformed into a PDM 
one, given in equation (|2.4p . via a change of variable 



u(a; x) = av (a; x) + 6, v(a;x) = J \J M(a; x') dx' , (3-2) 

and a change of function 

<j> n (u(a;x)) oc [M(a; x)]' 1/4 ip^ (x). (3.3) 



In (13. 2p . a and b are assumed real. The effective potential, defined on a possibly different 
interval, and the energy eigenvalues of the PDM Schrodinger equation are given in terms of the 
potential and the energy eigenvalues of the constant-mass one by 

M" 7M' 2 

Vett (a; x) = a 2 U(av(a;x) +b) + -^-^ - + c (3.4) 

and 

= a 2 e n + c, (3.5) 

where, as before, a prime denotes derivative with respect to x and we have introduced an additive 
real constant c, not present in [35] . 

We can reformulate the PCT, defined in (|3.2p and (13. 3D . in terms of the deforming function 
/(a; x) as 

/" x fix' 

u(a; x) = av (a; x) + 6, v(a;x) = / -77 jr, (3-6) 

7 f{oc;x') 

0„(«(a;x)) oc (3-7) 
Then, on taking (|2.5p into account, equation (|3.4p is replaced by 

V(x) = a 2 f/(aw(o;;2;) + 6) + c. (3.8) 

In most applications, this PCT is used in the following way (sec. e.g., |36j). One starts from 
a given exactly solvable constant-mass Schrodinger equation, hence from some known t/(w), E n , 
and 4> n (u). One chooses a PDM and some parameters a, b, which means some change of va- 
riable (|3.2p . As a result, one obtains an exactly solvable PDM Schrodinger equation, containing 
an effective potential given by (|3.4p . and whose eigenvalues and eigenfunctions can be directly 
derived from equations (13. 5p and (13.31) . respectively. There are variants of this approach, wherein 
the dependence of the variable u on the PDM is different from that given in (|3.2p . but otherwise 
the method remains the same (for some examples see, e.g., |37|). 

Here, in contrast, we are going to make use of our acquaintance with M(a;x) (or f(a;x)), 
V e ff(a;x) (or V(x)), and Et\ coming from the DSI method, to determine tpn(x). From 
the dependence of on n, we first guess which type of potential U (u) in the constant- 

mass Schrodinger equation (|3.ip may give rise to such a dependence through equation (|3.5p . 
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Furthermore, from the latter equation, we determine the values of the constants a 2 and c, as well 
as all or only some parameters of £/(it)0 In a second step, we obtain the change of variable (|3.6|) 
from our knowledge of f(a;x) and \a\, the sign of a and the value of b being determined from 
the known domain of definition of U(u). The third step then consists in checking equation (|3.8p 
and calculating the remaining unknown parameters of U(u), whenever there are some. Finally, 
since the wavefunctions 4> n ( u ) of equation (|3.ip are explicitly known, we deduce tp^*\x) from 
equation (|3.7p . 

Before listing the results for ^^(x), we shall proceed to illustrate our method with two 
examples: a simple one (corresponding to the Poschl-Teller potential), wherein the PCT does 
not change the nature of the potential, and a more complicated one (corresponding to the 
Morse potential), wherein the PCT deeply modifies the potential and the number of bound 
states acquires an intricate dependence on the parameters. 



3.2 Poschl Teller potential 

From Tables [1] and [3l we observe that for the deforming function f(a;x) = 1 + asm 2 x (— 1 < 
a / 0), the Poschl-Teller potential V(x) = A(A — 1) sec 2 x, defined on (— 5 , has a quadratic 
bound-state energy spectrum made of an infinite number of levels, 

£l(«) = (A + nf - a (A - n 2 ) , n = 0,1,2,.... (3.9) 
On the other hand, it is well known [38] that the same type of potential 

U(u) = A 1 (A! - 1) sec 2 u, (3.10) 
defined on (— ■§, §) and used with a constant mass, has a spectrum 

e n = {A' + n) 2 , n = 0,1,2,..., (3.11) 

with essentially the same characteristics. Equations (j3.5[) . (|3.9p . and (|3.1ip then lead to the 
relations 

a 2 = l + a, A' = -^—, c = — — A(A - 1 - a) = — —A(A - 1), (3.12) 
1 + a 1 + a 1 + a 

where in the last step we used the expression of A given in Table [2[ 
The second equation in (|3.6p now yieldH 

v(a;x) = arctan (y/l + atanx) . 

V 1 + a 

Furthermore, it is obvious that if we assume a = \f \ + a and b = in the first one, the variable 



u(a; x) = arctan (\/T+~a tan x) (3.13) 
is defined on (— §, f), as it should be. From the inverse transformation, we get 

tan u ,„ „ , . 

tanx = . (3.14) 

1 + a y ' 



1 It may happen for some potentials U(u) that not all their parameters appear in their eigenvalues e„. For the 
Scarf I potential, for instance, e n only depends on A, but not on B. In such cases, equation (|3.5|l is not enough 
to entirely determine U(u). 

2 Here and in other cases, we use the principal value of all inverse trigonometric functions. 
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On inserting f|3. 14j) in the definition of the starting Poschl-Teller potential, we arrive at the 
expression 

V(x) = — — (sec 2 u + a) , 
1 + a 

which agrees with equation (|3.8p when we take equations (|3.1Up and (|3.12p into account. 

Since, for constant mass, the bound-state Poschl-Teller wavefunctions can be expressed in 
terms of Gegenbauer polynomials as 



4> n {u) oc (cos u) A ^(sinu), n = 0, 1, 2, . . . , 
the use of equations (|3.T|) and f|3. 13j) leads to the searched for wavefunctions 

^\x) oc [/(a; x))-^ A ' +1 \^x) A 'C^\t), n = 0, 1, 2, . . . , 

/ 1 + a . . .n 

t= \ T? \ ' f{a;x) = 1 + asm x, 

V f{a;x) 

for the Poschl-Teller potential and the present PDM. Here A' should be replaced by its expression 
in terms of A and a, given in equation (|3.12p and in Table [5J 

As a final check, we may compare the obtained wavefunctions with those coming from the 
DSI method and given in Tables 2] and Knowing the former indeed provides us with a hint for 
solving the differential-difference equation satisfied by the polynomials P n (X;y). The changes 
of variable 



1 + Q 



l + (l+a)y 2 ^(l + a)(l-i 2 ) 

and of function 

P n (X;y) = ^ A 'Hi-tT n/2 c( A '\t), (3.15) 

(A 1 ) 

where 7n is some constant, transform the equation fulfilled by P n into the backward shift 
operator relation for Gegenbauer polynomials [39], 

- '"T'l'"" ^" 1 " = ( (1 " e) s + (1 " 2a)t ) cf){t) - {3M) 

where a is replaced by A' . In (|3.15p . we have to choose 

7 y> = q+ a y /2 +w*o 



2 ri r(2^ / + n)r(A' + n)' 
This completes the derivation of wavefunctions for the Poschl-Teller potential. 

3.3 Morse potential 

From Tables[H El andO we observe that for the Morse potential V(x) = B 2 e~ 2x - B(2A + l)e~ x 
(A, B > 0) and the deforming function f(a\x) = 1 + ae~ x (a > 0), defined on the real line 
-co < x < oo, the bound-state energy spectrum is made of a finite number of levels, whose 
energies can be written as 

E w = -„; = --LA- g!g±^±l)l V „ = ,i,. ..,„_, ,3.17) 
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where A„ = A — na and n max is such that 



"max 

(2|A| + n max a) < 2\X\fi < (n max + 1)[2|A| + (n max + l)a\. 
Such an energy spectrum looks like that of an Eckart potential 

U(u) = A' (A' - 1) csch 2 u - 2B' coth u, < u < oo, 
for a constant mass, which is given by |17j 

B' 



(3.18) 



-{A' + nf 



n 



A' + n / 

i A ' + ? Vax) 2 <B' < (A' + n' max + l) 2 . 



1 ??' 



By comparing (13.170 and (13.190 with (13.50 . it is indeed straightforward to obtain 



1 

4 : 



A' 



|A| 



a 



1 



a- 



B[B + a(2A + l)], 



(3.19) 
(3.20) 

(3.21) 



Furthermore, with such parameters A' and B' , condition (|3.20p directly leads to restriction (|3.180 
if we set n' miiy . = n max . Hence, the latter, whose origin in [TB] was based on the behaviour of the 
wavefunctions ^^(x) for x — > ±oo, gets here a very simple derivation in terms of the PCT. 

The remaining steps of the procedure work as for the Poschl-Teller potential. We successively 
obtain 



i 

2' 



v(a; x) = ln(e x + a) 
leading to 



€ X + OL 

< u(a; x) = - In < oo, 

2 a 



p In a, 



and 



^)(x)oc/-k^ + " +i+ ^: 



cxp 



A' + n 



B' 



X Prl 



-A'-n+ 



B' 
A' + n 



-A'-n- 



A'+nj 



(t), 



A' + n, 
t = l + 2ae~ x , 



(3.22) 



where in the last step we used the known wavefunctions <j) n (u) of the Eckart potential for constant 
mass in terms of Jacobi polynomials |17j . In (|3,22p . A' and B' should be expressed in terms 
of A, B, and a, as shown in (13.210 and Table [2l 

In the present case, the differential-difference equation for P n (X,u;y), given in Table 
amounts to a combination of the recursion and differential relations for Jacobi polynomials (see 
equations (22.7.1) and (22.8.1) of @Q|), 



2(n + l)(n + a + /3 + lj_ p ( a ,p) ^ 



2n + a + (3 + 2 
d 



n+l 



(1 - t 2 )— + {n + a + P + l)[t + 
at 



a — (3 



2n + a + + 2 



(3.23) 



with q and [3 replaced by —A' — n — 1 + A ,+ n+1 and —A' — n — 1 



B' 



A'+n+l 

relation between both approaches is obtained through the transformation 

, ... (-A'-n+-M—-,-A'-n--^—) 



, respectively. The 
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(A / } _ n!r(2 J 4' + 2n)r(,4 / ) 
7n ~ 2™r(2,4 / + n)r(A' + n) ' 

It is worth observing that the expression (j3.22j) . obtained here for the bound-state wavefunc- 
tions ipn\x), differs from equation (3.16) of [16], which is the result of a PCT starting from 
the PDM radial oscillator Schrodinger equation. In our previous result, the Jacobi polynomi- 
als indeed depend on the variable (— 1 + ae~ x )/(l + ae~ x ) instead of t = 1 + 2ae~ x , while 
the other factors are also slightly modified. On writing the Jacobi polynomials in terms of 
hypergeometric functions and using a linear transformation formula for the latter (see equa- 
tions (15.4.6) and (15.3.4) of [40]), it is however straightforward to check that both results are 
in fact equivalent, as it should be. 



3.4 Results 

For the potentials and deforming functions of Table H] we list the changes of variable (|3.6p . 
the potentials U(u) used in (|3.ip . and the resulting PDM wavefunctions ipn (x) in Tables [6] [7] 
and El respectively. 

Table 6. Variables u(a; x) = av(a; x) + b and their domain of definition. 



Type v(a;x) 



Domain 



SHO a 

SHO fe 

RHO 

C 

M 

E a 
E 6 
PT 

S 

RM 



1 , ax+8 

— g arccot — j^- 
-^arccot (*-£)] 

arctan(y / ax) 
±ln(l + ax) 
ln(e x + a) 
I arccoth 



coth x+ ^ 
-1) 



£ <%L CCU III J — 

(coth x — 1) _1 
,}, arctan (v/1 + a tanx) 

yfl+Oi \ ' 



fl+a 

vT-q 
j arccot 



tan §+a 

cot 



arctan ta , n 2 = 

V 1— or 



-5 







7T 

2 


2^ 




Q 

2 





1 




2 









1 





v 7 ! + " 





Vl-a 2 




S 






1+a 



< U < 7T 

-§<«<! 

-§<«<§ 

< u < oo 
< u < oo 

< u < oo 
< u < oo 

-§<«<§ 

-f <u<f 

< U < 7T 



In Table [8] we remark that in addition to Laguerre, Gegenbauer, and Jacobi polynomials, 
there appear the less known Romanovski polynomials R^*'^\t) [H] for the shifted harmonic 
oscillator and the Rosen-Morse I potentials. We indeed follow here a recent approach to the 
wavefunctions of the latter potential for constant mass [32], which has stressed the interest of 
employing such polynomials, solutions of the second-order differential equation 

((1 + t 2 ) ^ + (2{3t + a)j t - n(n - 1 + 2/3)) Rfr®(t) = 0, 

instead of Jacobi polynomials with complex indices and complex arguments, as is usual [17] . 
For more details on the properties of Romanovski polynomials, we refer the reader to [43] . 

Finally, in Table EE we provide the connections between the wavefunctions listed in Tables H] 
and [H To establish those results, we have used not only equations (|3.16p and (|3.23j) . but also 
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Table 7. Potentials U (u) and additive constants c. 



Type U(u) 



SHO a A(A - 1) esc 2 u + 2B cot it 
x 



(a- 25 2 )^ + ^dio + d 



A = -, B 



r! I Pu 
2a8 \ 2a 



SH0 6 A(A - 1) sec 2 u 



-Aad 2 



Act 



A = - 

a 



RHO [A(A - 1) + B 2 ] sec 2 u-B(2A-l) sec u tan u 

A = J(L + l + g),fl = J(i+l-S) 
^(^ - 1) csch 2 u-2B coth it 

A = L + 1, B = ^- + L(L + 1) 
A' (A' - 1) csch 2 u - 25' coth it 
A' = -A B' = £fl[J3 + a(2A + l)] 
A(t4 - 1) csch 2 u - 2B' coth it 
B' = \(B + \aA(A - 1)) 



a L(L + 1) + 



C 



M 



E 



E 6 



2Z , 



la- 

f [2Z + aL(L + 1)] 
-^B[B + a (2A + l)] 
a[B + 5A(A- 1)] 
-2B 



aA(A-l) 
1+a 



U M z 

Z = B-A(A-1), L = A-l 
PT A 1 (A 1 - 1) sec 2 u 

1+a 

S [A'(A'- 1) +5 /2 ]sec 2 n- J B / (2A' - 1) sect* tan u -j^{a[A(A - 1) + B 2 

A' = |(1 + Ai + A 2 ), B' = i(A 2 - Ai) + B(2A - 1)} 



RM A(vl - 1) esc 2 it + 2B' cot u 
B' = \[B- \aA(A - 1)] 



- /?) A(A - 1) - aB 



a combination 



(„ + !)(« + n)L^; 2 (t) = ((a - l)t| + a(a - 1) - (a + n)tj L^(t) 



(3.24) 



of the recursion and differential relations of Laguerre polynomials (see equations (8.971.5) 
and (8.971.3) of [44] ) . the backward shift operator relation for Jacobi polynomials (see equa- 
tion (1.8.7) of [39]), 

2(n + l)P£f(f) = (- (1 - t 2 ) ± + a - /? + (a + + 2)f) P^ +1 ^(t), (3.25) 
and the differential relation of Romanovski polynomials (see equation (41) of [43| ) . 



n + 2(3 - 1 (q.,/3) 
2(n + /3) 



fi+i 



1 + r) — + (n + 2/3 - 1) ( t + 



a 



2(n + p) 



R^\t). 



(3.26) 
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Table 8. Bound-state wavefunctions resulting from the PCT approach. 



Type Vi a) (x) 



SHO a 
SHO b 

RHO 

C 
M 

E a 

E 6 

PT 

S 

RM 



/ _i (A +n+l) exp 



1^- arccot a 



A+ 



r^ A+1) d A \t) 

xA+ B r A-lpiA-B-h^B-l) {t) 



A+n , 



(/) 



/ 2 V A +n y g x p 



x P„ v 



-A'-ri+ 



B' 
A 7 ^ 



-A'-n- 



B' \ 



B' 



A'+n 



(cothx + 1)1 (cothx + 1 + a) 2 ( A+n+1+ A+«, 



1 ( A-t-r, -^M (-A-n+-^--A-n--£— 
x (cothx - l)"H A+n -^Jp n V A+n ' A+n > ( t ) 



(n+L+l)(cothai-l) 



(cothx + 1)2 (cothx — 1) L 2 exp 
xLf +1) (i) 

/-|(Ai+A 2 +2) (l _ sinx) I(A 1 + I) (1 + gin;r) |(A 2 + I) 
xPi Al ' A2) (t) 



■§(A+n+l) 



(sinx) 



A+n 



cxp 



' arccot cot x+ 2 



1 A+n 



X R r \ 



A+n ' 



-A-n+1 



ax+f} 



(x 



— 1+ax 2 
1+ax' 2 

2+ax 



U! ) 



1 + 206-* 



2 coth a+q 
2<5 



2Z 



(n+L+l)(cotha;-l) 



1+q 



smx 



sm x+a 
l+a sin a; 



2 cot x+a 
28 



4 Conclusion 

In the present paper, we have completed the study of a class of exactly solvable PDM Schrodinger 
equations, undertaken in [15] and pursued in [161 127j . by constructing their bound-state wave- 
functions in general form. This has been made possible by combining the previously used DSI 
method with PCT connecting the PDM Schrodinger equations with constant-mass ones with 
similar spectra. 

It should be stressed that it is only the association of both approaches that has allowed us 
to build and to fully solve the former equations. As shown in Table [H the complexity of the 
changes of variable involved in the PCT would indeed make it rather unlikely to guess them 
in order to directly construct the PDM equations from the constant-mass ones. The previous 
determination of the PDM and of the corresponding spectra by DSI techniques has therefore 
played an essential role. 
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Table 9. Changes of variable and of function relating the wavefunctions arising from the DSI and PCT 
approaches and equations used in the comparison. 

Type t P n (X;y) Equation 

SHO a 
SH0 6 
RHO 

C 

M 

E a 

E 6 
PT 
S 

RM 
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